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1. INTRODUCTION 

The concept of an intelligent evacuation system involves information systems in affected 
placements in several alternative locations. This concept is a solution to dynamic evacuation activities that 
will be carried out in mitigation activities [1]. Several evacuation locations were competed to determine the 
strategy to be used as an evacuation location by considering the geographical aspect where the impact of the 
disaster could not be predicted. The purpose of determining the evacuation location is to maximize the 
number that can be saved and the disaster's effect may well be reduced. Several aspects are considered in the 
evacuation process, an example of this would be the amount that could be mobilized, how many people can 
be transported at a time, and utilized to get to an evacuation site by the quickest way possible [2], [3]. 

In managing evacuation, the area has a constraint function, such as transporting and concentrating 
those affected by the disaster to a safe area [4], [5]. The evacuation scale is an important part that is 
considered in evacuation activities. In general, evacuation consists of small-scale evacuation and long- 
distance evacuation. Small-scale evacuation is usually guided by emergency evacuation immediately and 
emergency conditions that are affected by a small space, for example, the occurrence of a small explosion, 
the collapse of a building, or the occurrence of a fire. Long-distance evacuation generally requires a vehicle 
as a mode of transportation and is generally instructed by disaster-related institutions or related ministries 
that communicate with other units [6], [7]. Evacuation activities are required to move quickly to place the 
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affected victims in the disaster location to a safe location by considering the absence of congestion and chaos 
in the traffic situation throughout the evacuation process [8]. 

The evacuation process requires an evacuation route that is regulated in the concept of transportation 
management. In the evacuation process, it is immediately possible to find the fastest route, such as the pattern 
of solving the traveling salesman problem (TSP) if the starting and ending nodes are the same and each node 
must be visited and is also called the vehicle routing problem (VRP) [9]. The selection of the fastest route can 
also be solved by using the open vehicle routing problem (OVRP) method. OVRP, on the other hand, does 
not require the vehicle to be returned to its original pool once it has visited the last customer on the route 
[10]. OVRP problem-solving uses heuristic or metaheuristic methods such as tabu search [11]. Previous 
research that has been done discusses the application of the disaster traffic management model with the 
heuristic open vehicle routing problem with time window (HOVRP) in a contra-flow state [12]-[14]. 

The disaster evacuation handle routing time model is an optimization of disaster evacuation 
transportation management conditions contra-flow not to occur in all route sections with the central 
assumption that safe areas have been determined with an optimum time of 400 minutes during the evacuation 
process [15]-[17]. In this study, we will discuss the dynamic model of determining the location of disaster 
evacuation using game theory and determining the effective evacuation route. Game theory applications can 
be used in various emergencies [18]. One of the game theory developments is the combinational game, which 
is looking for an optimal strategy in several possibilities [19]. 


2. PROBLEM STATEMENT 

One of the techniques used to produce optimum decisions is linear programming, a mathematical 
method that places limited resources for optimal outcomes or results, such as maximizing profit or 
minimizing costs [20]-[22]. In linear programming problems, the solution can be in the form of fractions to 
produce the optimal value in the form of real numbers. The result of the solution may deviate from the goal 
because the solution is rounded to the nearest integer. Not a few problems in real life require solving decision 
variables in integer form and problem-solving models must be sought to obtain optimum integer solutions. In 
integer programming, some or all of the decision variables are integers as part of linear programming 
development [23], [24]. Evacuation problems can be done by using a linear program approach, which is 
game theory. In its implementation, it can be found in activities that are competitive in nature which allows 
competition or conflict to occur. The impact of this competition occurs on two people or two parties or a 
number of groups of people or groups [19]. 

The game is a form of competition between 2 people or parties or between 2 groups or teams that 
are close together and use rules that are known by the adjacent parties [25], [26]. On the other hand, the 
philosophy of game theory is a mathematical approach that defines an atmosphere of competition and 
conflict between various needs [27]. Using this concept, we can better understand how various competitive 
environments influence decision-making and how two or more requirements are interconnected. 


2.1. Game theory on disaster management 

The strategy used in determining the evacuation location uses a mixed strategy. In this strategy, a 
player or company will use a mixture of more than one strategy to get maximum results, where: x; = 
Probability of player A with strategy i(i = 1,2,3, ...,m), yi = The probability that player B chooses strategy 
JG = 1,2,3, ...,n) [28], [29]. So that; 

i=1%j = Aj- = 1 (1) 

Xi, yj 2 0, for every i and j (2) 
The minimax criteria for mixed strategy are as shown in: player A has: 

xi (xi 2 0, Xx = 1) (3) 
which will produce: 

max x; {min (Yih 1 1Xi, Bi=1 li2Xi s => Lier AinXi)} (4) 
Player B has; 

yO = 0, Ray = 1) (5) 
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which will produce: 
min y; {max (i; Cj 02.721 02jYj s => Deja Anji) f (6) 
EV += Vita Dj- ayXi Yj (7) 


The payoff value obtained from the subtraction between player 1 and player 2 is assumed to have a value as 
in the matrix below. The strategies used are as shown in: 


EL, =Loc 1 

EL, = Loc 2 

EL, = Loc 3 

EL, =Locn 

Xa = Place elevation 

Xp = Total population 

Xe = Evacuation route length 
Xa = Hierarchy of the roads 
Xe = Type of transportation 
Xf = Evacuation area 


Using players and tactics, the optimum solution to the return calculation is found. In deciding where 
to evacuate, the tactics used throughout the evacuation may be advantageous. The computation of each player 
is expressed by associating Loc 1 (EL,) with Loc 3 (EL3), Loc 2 (EL) with Loc 3 (EL3), and Loc 1 (EL,) with 
Loc 2 (EL,). At each evacuation location, a strategy can be obtained that can be used to achieve an advantage 
in the evacuation, which is assumed to have the value of the game player comparison as follows. It is assumed 
that the optimum solution is the optimum solution for the Loc 1 and Loc 3 matrix as shown in Table 1. 


Table 1. Optimum solution of location 1 and location 3 matrix at location 1 
X, Xis X- Xa Xe Xp Z 
o 001 0 0 0 0.03 0.04 


So: 


1 


V = — = 
0,04 


25 


Xa =Xa*V=0*25=0 


X, = Xp * V = 0.01 * 25 = 0.25 


X. =X, *V=0*25=0 
Xa =Xq*V =0*25=0 
Xe =X,*V =0*25=0 


Xp = Xp * V = 0.03 * 25 = 0.75 


Because K=20 has been increased by the acquisition matrix components shown above, the game's 
worth rises to V = 25 — 20 = 5. Asa result, the best course of action is 5. The best course of action is obtained 
for Loc 1, namely the strategy for population and evacuation area, with a game value of 5 as shown in Table 2. 


Table 2. Optimum solution for location 1 and location 3 at location 3 
Ya Yw Y Ya Yo ý Z 
0.01 002 0 0 0 001 0.04 
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So: 


Y, = Y, * V = 0.01 * 25 = 0.25 


Y, = Y, * V = 0.01 * 25 = 0.25 


Y,=Y¥,*V =0*25=0 
Y4 = Y4 * V = 0.01 * 25 = 0.25 


Y, = Y, * V = 0.01 * 25 = 0.25 


Y, =Y *V =0*25=0 


The game value is now V = 25 — 20 = 5 because of the addition of the acquisition matrix 
components to K=20, as a result, 5 is the optimal solution. The best course of action is achieved for Loc 3, 
that is the altitude of place, total population, hierarchy of the roads, and type of transportation strategies with 
game value 5. It is assumed that the optimum solution is the optimum solution for the Loc 1 and Loc 2 matrix 
as shown in Table 3. 


Table 3. Optimum solution of location 1 and location 2 matrix 
Y, Y, Minimum 
Xa 7. 3 3 
X; 7 83 7 
Maximum 7 13 


The value of 7 is seen in the table above, if a saddle point exists, then the game is in its most ideal 
state. The optimal strategy for location 2 is the population size strategy and the optimal strategy for Loc 1 is 
the evacuation area strategy. It is assumed that the optimum solution is the optimum solution for the Loc 2 
and Loc 3 matrix as shown in Table 4. 


Table 4. Optimum solution of location 2 and location 3 matrix at location 2 
Zo Xo Xo Xa X- X Z 
0.01 0 0 0 001 0.02 0.04 


so: 
1 

V= 004 =25 

Xa = Xa * V = 0.01 *25 = 0.25 

X, =X, *V =0*25=0 

X, =X,*V =0*25=0 

Xa =Xq*V =0*25=0 

Xe = Xe * V = 0.01 *25 = 0.25 

Xp = Xp *V = 0.03 * 25 = 0.75 

Because of the addition of the acquisition matrix components K=40, the game's value is 
V = 25 — 40 = —15. As a result, -15 is the best course of action. The best course of action for Loc 2 is 
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achieved, that is the strategy of place height, type of transportation and area of evacuation place with a game 
value of -15 as shown in Table 5. 


Table 5. Optimum solution of location 2 and location 3 matrix at location 3 
Xa Xp © Xa Xe X Z 
0.03 0 0 o 001 0 0.04 


so: 


1 


Y= 04 = 


25 


Y, = Y, * V = 0.03 * 25 = 0.75 
Y, =Y¥y*V=0*25=0 
Y.=Y,*V=0+*25=0 
Ya =Y¥gxV =0*25=0 
Y, = Y, * V = 0.01 * 25 = 0.25 
Y, =Y;*V =0*25=0 


The game value is now V = 25 — 40 = —15 because of the addition of the acquisition matrix 
entries to K=40, as a result -15 is the optimal solution. The best course of action for Loc 3 is achieved such as 
altitude of place and type strategy of transport with a game value of -15. Accomplishing the goals of the 
evacuation procedure is dependent on the evacuation site. In determining the evacuation location, one must 
provide alternative locations due to the impact of unpredictable damage. Game theory, a combinational game 
to find the optimal strategy in several possibilities, can be used to determine alternative evacuation locations. 
The determination of this evacuation location is given a set of N, which will then be optimised in selecting 
an evacuation route as a set of final spots or safe destinations N, = {vn we Ung Vas} In solving the problem 
of choosing an evacuation location, it is assumed that it has more than 1 location or is given a set of L, with 
influencing variables such as altitude (X,), population (X,), length of evacuation route (X,), road hierarchy 
(Xa), type of transportation (Xe), and the area of the evacuation site (X-). The simulation of determining the 
evacuation location can be seen in the following Figure 1. 


Figure 1. Simulation of alternative evacuation locations 
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2.2. Optimisation of evacuation routes 

After the evacuation location is determined, the evacuation route will be determined. Evacuation 
route selection is part of transportation management. In determining the evacuation route, a mathematical 
model is made to minimize the time in the evacuation process. The flow of vehicular traffic at the time of a 
natural disaster can be expressed in a directed graph G = (V,E) with V = {v%,¥4,...,Un} is the set 
of evacuation point spots. Node vo is the starting point for evacuation, with a record of the number of 
nodes vg get more than one. It is assumed that at these M nodes, the vehicle fleet is located. Then 
E= {(v, vj): vp vj EV LF TJ + 0} is the set of vehicle paths. In the VRP, all vehicles must depart from the 
central node and return to that node. Whereas in OVRP the vehicle does not return to the central node. This is 
reasonable because the flow of the route starts from the evacuation center Vo and goes to the safe area node vy. 

The model that will be generated follows the VRP model pattern. The type of vehicle used for 
evacuation is heterogeneous. So, the current model pattern is heterogeneous VRP. For evacuation it is 
necessary to have time windows. Therefore, in the realm of VRP, this condition refers to the heuristic open 
vehicle routing problem with time window (HOVRPTW) model pattern. 

A route is used to symbolize every arc in E (i,j), linking nodes i and j through a route. This kind of 
network is known as a static network since every arc in the network only connects part of the nodes in the 
network to each other. A parameter is attached to each route and node. As a network, each node k represents a 
place with an initial population pẹ and capacity qx. Given a capacity of c;j for every path (i, j), where (i, j) € 
E. A path's capacity is the amount of flows it can handle in a given period of time, provided there is no 
congestion. It is the total of vehicles per hour on each lane that determines capacity on a lane-based roadway 
network. For each trajectory, it defined travel time t;; where (i, j) E€ E. If the road is clear of evacuation, it is 
considered that the t;; constant represents the average speed of travelling over the arc (i,j). It is commonly 
known to as the free flow rate or grace time (lead) for the bow by this particular characteristic (i, j). 

When it comes to how many people can be evacuated at a given moment, the line capacity is always 
considered to remain constant. Unfortunately, in reality, the capacity of the track is not always the same. It is 
the total of entities present in an arc that determines the capacity of that arc at any moment. Flow-dependent 
capacity reduces the number of network flow issues, and adds new limitations, by incorporating it into the 
system. The objective to be achieved from this problem model is planning to regulate the flow of vehicles 
from the initial location of the disaster to a safe location so that the vehicle's travel time is minimal. The 
notation used is as shown in: 


Set: 

N : Set of nodes 

K : Vehicle set 

No : Evacuation starting node set, No = {Vo1, Voz, +) Voa} 

N, : Transfer node set N\ {vo}, {Vn} 

N, : Set of end nodes (safe destination), No = {Vn1; Vnz =» Uns} 
T : Travel time set 

Parameter: 

Qk : Vehicle capacity k € K 

tij : Travel time from node i € Np to node j € N, 

[aibi] : The earliest and the last time on the node i € N, 

Cij : The path capacity of the node i € Nọ U N; to node j E€ N, UN, 
St; : Service time at node i E Nọ UN, 

td; : Slowest time to arrive at node i E€ N, 

3. METHOD 


The research was conducted using descriptive research methods. This research was conducted to 
determine the value of the independent variable, both 1 or more independent variables and no comparison 
was made between several variables [15]. The data collection technique used in this research is a literature 
study that collects data by studying various literature books and related documents to find solutions to 
problems in determining several alternative evacuation locations and optimal routes or routes for evacuation. 
The model development carried out is a dynamic model determining the location of evacuation in disaster 
mitigation which is generally seen in the following Figure 2. 

The research starts from formulating a description of the problem, namely, how to determine an 
evacuation location based on the geographical conditions of an area. Next, describe the model for 
determining the evacuation location and determine the variables and parameters of the problem. Next, the 
game theory model will be used to solve the issue of locating the evacuation site, and once the findings are 
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obtained, the development of the model for locating is discussed. The next stage is the determination of 
model variables and parameters to complete the evacuation route or transportation route and followed by 
making an optimization model of the evacuation route/path and validation by testing the data using linear, 
interactive, and discrete optimizer (LINDO) software and after being valid it becomes a new model for 
determining evacuation locations in disaster mitigation. 

- Problem parameters 


Formulation of 
Problem Description 


- Variable Description of the Game Theory Model for 
- Model Establishment of Determining Evacuation 
Parameters Evacuation Site Locations 
Determination Model 


Evacuation Route m 
Optimization Model m Model Validation = New Model | 


Figure 2. Research methodology 


Description 
Determination of 
evacuation location 


- Variable 


4. RESULTS AND DISCUSSION 


Intermediate nodes signify where evacuation flows are converging (merging) or crossing. 
Minimizing the amount of time needed to evacuate in order to salvage as much material as possible is the 
primary objective. The objective function is expressed in the following expression. 


Min = Wien, È jen, tij Lex xf + Vien, Lien, i+; tij Leek x + Pieno Ljen, tij Lex yis + Vien, Lien, i+; tij Leek yÉ (8) 


4.1. Mathematical model formula 
With obstacles or requirements that need to be met: 


È kek Licno xi = d,Vj EN, UN? (9) 


This equation ensures that the vehicle for each type departs from the l node of the initial disaster location to 
the next j node in the E path. 


jen, x5 =K Vrek viuN, (10) 


This equation is to ensure that the number of vehicle trips departing from the initial disaster location 
is at most K paths. 


Baene Lex Vai + Lyen,;i#j Leen Xf =1 Wie Ny (11) 


This obstacle requires that a path cannot be traversed unless it starts from the initial location of the 
disaster. 


Djen, Xi? = Djen, Xi; Wk E€ K, Vi € No, Vt ET (12) 


This expression is to force that the vehicle will visit and leave a node at the time of travel t. 
Ljun xg =1 Vk EK, Vi € No (13) 


This equation is to ensure that for each vehicle there is only one path from the selected initial 
disaster location. 
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Dien xika < L Yj EN UN, Wk EK, VteT 
There are artificial nodes (artificial). 

Erer Lien Xani = O Yj E Ny, Vt ET 
It is not allowed to start service from artificial nodes. 


tk i 
È kek L jeNiun xij = 0 vi EN, YtET 


(14) 


(15) 


(16) 


Vehicles that have left the initial location of the disaster are not allowed to return to the initial 


location. 
Ea pen,izj XE < ST-1, YkEK,YytET 

Sub-tour elimination process (ST). 
Lupeni<j Xj — Lapenj<ixi =0 YkEK 

This equation is to maintain the flow of vehicle flow to continue. 
Erer Lien Xi = 0, 

To prevent loops. 


k 
Vij S CijX 


ij Wie No, Vj EN, UN2,i #j,VK EK 


The equation to state that the number of path currents does not exceed the path capacity. 


IF Sa; Dien, x, Vie N Vk EK 
A; Dien, Xj < lF + ti < bi Djen, x5, Vi € No, Wk E€ K 
It should be noted the time of arrival at a node and the time of departure. 
xi vii € {0,1}  Y(Lj)E€N,Yk E€ K 
l¥>0 vieN,ykEeK 
Requirements for decision variables. 


4.2. Model trial 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


Trial of optimization models to determine evacuation routes or paths using LINDO software with 
the input means of transportation used is heterogeneous. It can function as a means of transportation that 
moves in the same direction and does not return to the starting point/node as the starting point of the moving 
vehicle with the following input. Based on the calculations, the minimum travel time of vehicle in matix X 


and Y are shown in Table 6 and Table 7, respectively. 


Table 6. Minimum vehicle travel time matrix X Table 7. Matrix of minimum vehicle travel time Y 


Ji XO XI X2 X3 X4 X5 Ji yO yl y2 y3 y4 y5 
i=1 5 0 0 0 0 0 i=l 5 0 0 0 0 O 
i=2 8 9 0 0 0 0 i=2 8 9 0 0 0 ODO 
i=3 4 10 6 0 0 0 i=3 4 10 6 0 0 0 
i=4 7 8 9 12 0 0 i=4 7 8 9 12 0 0 
i=5 6 12 11 10 8 0 i=5 6 12 11 10 8 O 
i=6_ 4 6 8 9. r-T i=6_ 4 6 8 9 12 7 
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Minimize 
5 X011 + 8 X021 + 4 X031 + 7 X041 + 6 X051 +4 X061 + 9 X121 + 10 X131 + 8 X141 + 12 X151 +6 
X161 + 6 X231 + 9 X241 + 11 X251 + 8 X261 + 12 X341 + 10 X351 + 9 X361 + 8 X451 + 12 X461 +7 
X561 +... 


Minimize 

5 Y011 +8 Y021 +4 Y031 +7 Y041 +6 Y051 +4 Y061 +9 Y121 +10 Y131 +8 Y141 +12 Y151+6 
Y161 +6 Y231 +9 Y241 + 11 Y251 +8 Y261 + 12 Y341 + 10 Y351 +9 Y361 +8 Y451 +12 Y461 +7 
Y561 


With the constraint function is as shown in: subject to 
X011 + X211 + X311 + X411 + X511 + X611 =1 
X022 + X122 + X322 + X422 + X522 + X622 = 1 
X033 + X133 + X233 + X433 + X533 + X633 = 1 
X044 + X144 + X244 + X344 + X344 + X644 = 1 
X055 + X155 + X255 + X355 + X455 + X655 = 1 

The results of the input above are as shown in: 
LP OPTIMUM FOUND AT STEP 23 
OBJECTIVE FUNCTION VALUE 
— 50.00000 

From LINDO’s calculation results, the maximum time to evacuate is 50 minutes with variable 
values and costs as shown in Table 8. In slack or surplus, it can be observed whether the resistance is active 
or not. If slack or surplus equal to 0 it means the constraint is active. Moreover, if slack or surplus not equal 
to 0 means the constraint is not active. For example, row 12 is an active constraint whose dual price is 
negative (-) i.e., 4, 000000 as can be seen in Table 9. This figure indicates that the accumulation of each part 
of the right-hand side of the obstacle will decrease the value of 4,000000. The output of the sensitivity 
analysis is as shown in Table 10. No. Iterations= 23. Ranges in which the basis is unchanged: obj. coefficient 
ranges. 


Table 8. Value of optimal decision variables 
Variable Value Reduced Cost 
X011 0.000000 1.000000 
X021 0.000000 4.000000 
X031 1.000000 0.000000 
X041 0.000000 3.000000 
X051 0.000000 2.000000 
X061 0.000000 0.000000 


Table 9. Slack or surplus 
Row __ Slack Or Surplus _Dual Prices 


2) 0.000000 0.000000 
3) 0.000000 0.000000 
4) 0.000000 0.000000 
5) 0.000000 -2.000000 
6) 0.000000 0.000000 
7) 6.000000 0.000000 
8) 6.000000 0.000000 
9) 6.000000 0.000000 
10) 6.000000 0.000000 
11) 6.000000 0.000000 
12) 0.000000 -4.000000 
13) 0.000000 -6.000000 
14) 0.000000 0.000000 


Table 10. Sensitivity analysis 
Variable Current Coef Allowable Increase Allowable Decrease 


X011 5.000000 INFINITY 1.000000 
X021 8.000000 INFINITY 4.000000 
X031 4.000000 0.000000 INFINITY 
X041 7.000000 INFINITY 3.000000 
X051 6.000000 INFINITY 2.000000 
X061 4.000000 INFINITY 0.000000 
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Based on the sensitivity analysis results above, it is explained that the variable column is the 
decision variable, and the current column current coef is the coefficient of the objective function. In the 
column, allowable increase on the variable X011 is the limit of the increase in value so as not to change the 
optimum value of the decision variable. In the allowable decrease column, the value is unlimited, which is 
the limit for the decrease, so that it does not affect the optimum value of the decision variable. 


5. CONCLUSION 

This research produces a dynamic model of determining evacuation locations in disaster mitigation, 
integrated with optimization of transportation routes related to the intelligent evacuation system. To 
determine the evacuation location with game theory algorithm and the determination of evacuation routes 
using the OVRP model can produce a new open dynamic model. The simulation model is tested using 
LINDO, which can minimize the evacuation time of 50 minutes from the evacuation route to a safe 
destination location with the assumption that there is no contra-flow to increase the population that can be 
evacuated. 
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